On the co-existence of matter dominated and accelerating solutions in f{Q)— gravity 
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Working within the theory of modified Gauss-Bonnet gravity, we show that FLRW-like power- 
law solutions only exist for a very special class of f(Q) theories. Furthermore, we point out that 
any transition from decelerated to accelerated expansion must pass through Q = 0, and no function 
f(Q) that is differentiable at this point can admit both a decelerating power-law solution and any 
accelerating solution. This strongly constrains the cosmological viability of /(£?)-gravity, since 
it may not be possible to obtain an expansion history of the universe which is compatible with 
observations. We explain why the same issue does not occur in /(i?)-gravity and discuss possible 
caveats for the case of /((J)-gravity. 
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INTRODUCTION 

One of the most remarkable developments in cosmol- 
ogy over the past decade has been the confrontation of 
the standard model of cosmology with observations, lead- 
ing to the conclusion that the Universe is accelerating in 
the current epoch. This rather counter-intuitive result 
has led to one of the most challenging theoretical puz- 
zles in 21st century physics, namely to provide an expla- 
nation for this late-time cosmic acceleration. Currently 
the most popular idea is that the energy density of the 
Universe is at the present time dominated by some mys- 
terious dark component known as dark energy and there 
have been many proposals offered to explain its nature. 
At the moment, the one which appears to fit all available 
observations (Supernovae la [H, Cosmic Microwave Back- 
ground Anisotropies [2| , Large Scale Structure formation 
Q, baryon oscillations QJ, weak lensing Q), turns out 
to be the so called Concordance Model in which a tiny 
cosmological constant is present [6[ and ordinary matter 
is dominated by Cold Dark Matter (CDM). However, de- 
spite its success, the A - CDM model is affected by signif- 
icant fine-tuning problems related to the vacuum energy 
scale, so other exotic negative-pressure fluids, often de- 
scribed in terms of scalar fields Q, have been proposed 
to address these issues. The problem remains however, 
that at the present time, there is no direct experimental 
evidence for the existence of the scalar fields responsible 
for the late time (and also the early time) accelerated ex- 
pansion rate of the Universe and this has led to a search 
for other viable theoretical schemes, many of which are 
based on the idea that the "dark sector" originates from 
modifications of the gravitational interaction itself. 

Currently, one of the most popular alternatives to 
the Concordance Model is based on modifications of the 
Einstein-Hilbcrt action. Such models first became pop- 
ular in the 1980's because it was shown that they natu- 
rally admit a phase of accelerated expansion which could 



be associated with an early universe inflationary phase 
Q. The fact that the phenomenology of Dark Energy 
requires the presence of a similar phase (although only a 
late time one) has recently revived interest in these mod- 
els. In particular, the idea that Dark Energy may have a 
geometrical origin, i.e., that there is a connection between 
Dark Energy and a non-standard behavior of gravitation 
on cosmological scales is currently a very active area of 
research. Additionally, such modifications are appealing 
due to the fact that they can evade constraints on the 
strength of the gravitational field which are restricted to 
below solar system scales. 

One of the most promising modifications to date are 
those based on gravitational actions which are non- linear 
in the Ricci curvature R and/or contain terms involving 
combinations of derivatives of R UM, EH, El- These 
theories have provided a number of ver 
suits on both cosmolog ical [3Eir 
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physical ID, 18 1 scales. One of the nice features of these 
theories is that the field equations can be recast in a way 
that the higher order corrections are written as an en- 
ergy-momentum tensor of geometrical origin describing 
an "effective" source term on the rig ht hand side of the 
standard Einstein field equations U 3] ■ In this Curva- 
ture Quintessence scenario, the cosmic acceleration can 
be shown to result from such a new geometrical contribu- 
tion to the cosmic energy density budget, due to higher 
order corrections of the Hilbcrt-Einstein Lagrangian. 

More recently modifications based on a Lagrangian 
density which is some general function of the Ricci 
scalar and the Gauss-Bonnet term f(R, Q) [H], [52] have 
been studied in the context of cosmology. In particu- 
lar a simple extension of Einstein gravity: f{R,Q) — 
R/2 + f (G ) has been investigated by a number of authors 
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231 [2J, |25|, |26j, |27[ . Corrections of this type can be mo- 
tivated from low-energy string effective actions and com- 
pactification of other higher dimensional theories. Uddin 
et al. [iH considered the existence and stability of power- 
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law scaling solutions in f(G) models and found that scal- 
ing solutions exist in the model f(G) — ±2y/aQ, where 
a is an arbitrary constant. Inspired by the work on /_( 
gravity by Carloni et. al. and Amendola et. al. 



Zhou et. al. [3l| used dynamical systems techniques to 



study the cosmology of f(G) models and argued that one 
could find cosmologically viable trajectories that mimic 
the ACDM cosmic history in the radiation and matter 
dominated periods, but also have a distinctive signature 
at late times 

In this paper, we investigate the cosmological viabil- 
ity of f(G) models by investigating the conditions under 
which one can find power-law solutions that mimic the 
standard FLRW expansion history of the Universe. We 
discover that such solutions only exist for a very special 
class of f(G) theories. Furthermore, we show that for 
there to be a transition between decelerated and accel- 
erated expansion phases, G must pass through zero, and 
there are no functions f(G) that are differentiable at this 
point that admit both a power-law decelerating solution 
and any accelerating solution. This seriously constrains 
the cosmological viability of /((?)-gravity, since it may 
not be possible to find an expansion history of the Uni- 
verse compatible with observations. 



FIELD EQUATIONS FOR GENERAL f(Q) 

We consider the following action within the context 
of four dimensional homogeneous and isotropic space- 
times i.e., the Friedmann-Lemaitre-Robertson- Walker 
(FLRW) universes with no spatial curvature, 



A 



-g[R + f(G) + C m ] 



(1) 



where R is the Ricci scalar, and / is general differentiable 
(at least C 2 ) function of the Gauss-Bonnet term, 



G = R 2 - AR a pR a0 



(2) 



and C rn represents the matter contribution. As we know, 
in four dimensions, the Gauss-Bonnet term is a total dif- 
ferential, and hence for f(G) = G, the field equations re- 
main invariant. However for other functions this term has 
non-trivial contributions to the field equations. For the 
homogeneous and isotropic spacetimes, the field equa- 
tions are: 

• The Raychaudhuri equation 



e 
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-y(p + 3P) + -0 3 /ssS 
(f-Gfg)-^Gfg-% 2 fg, 



(3) 



e - "" 3 

where is the volume expansion which defines a 
scale factor a(t) along the fluid flow lines via the 
standard relation = 3a/a, and f n g abbreviates 
d n f/(dG) n for n = 1,2, 



• The Friedmann equation 
8 
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^6 3 fggG + 3f = 3K 2 p + 3Gfg, 



(4) 



The total trace of the Einstein equations 
- 40 2 - 60 = 3k 2 (3P - p) + 12 (/ - Gfg) 



+ 36 ffggG 
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£/s + 80 2 /s, (5) 



• The conservation equation for standard matter 

p=-Q( P + P). (6) 

For FLRW spacetimes, the Gauss-Bonnet term is given 
by 
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(7) 



We can easily see that accelerating models have G > 0, 
while decelerating models have G < 0. In particular, any 
expansion history evolving from deceleration to acceler- 
ation must pass through G = 0. This observation will be 
of importance in the following discussion. 



Requirements for the existence of power— law 
solutions 



Let us now assume there exists an exact power-law so- 
lution to the field equations, i.e., the scale factor behaves 



a(t) = a t n 



(8) 



From now on, we assume that m is a fixed real number. 
If < m < 1, then the required power-law solution is 
decelerating, while for m > 1 it is accelerating. Since we 
know that within the standard paradigm, the expansion 
history of the universe underwent a power-law deceler- 
ating phase, it is important to study these kinds of exact 
solutions in our modified gravity models. 

We further assume that matter can be described by 
a barotropic perfect fluid such that P — wp with w G 
[—1,1]. From the energy conservation equation, we ob- 
tain 



p(t) = p t- 3m ( 1+w \ 



and the Gauss-Bonnet term becomes 



G = 24m 3 (rn - l)i~ 4 = a m t~ 



(9) 



(10) 



The negative sign of G for all decelerating models is 
reflected by a m < for the power-law models with 
< m < 1. 
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Using the background solutions above, we can write 
the Friedmann, Raychaudhuri and trace equations in 
terms of functions of time t only. We can assume with 
no loss of generality that t > 0. We then solve (flT)]) for 
t and re-write these equations in terms of G, f(G) and 
its derivatives. The Friedmann equation for example be- 
comes 

,- M+3m ,JT- «f«w - K /_e_\S"*< 1+ »> _ o 

V am a m \a m J 

(11) 

where K = pq/o^ 1+w \ Note that for the power-law 
solution |8]). G/ot-m is positive at all times by definition 
(|10p. and therefore equation (|lip is real- valued over the 
range of G- 

Since we want © to be a solution at all times, i.e. 
G spans over an entire branch of the real axis, we can 
interpret (jTTJ) as a differential equation for the function 
/ in G space. This differential equation has the general 
solution 

rz _ 3 , , 1 m 

f(G) = A m \jG+B mw Gi m{1+w) +C 1 G+C 2 Gi-~ . (12) 

We have abbreviated G = 24Q/a m , which we repeat is 
positive for the power-law solution l[5]). and 
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B n 



m(l-\-w) 



(to- l)K 



4 - to(19 + Ihw) + 3to 2 (4 + 7w + 3w 2 ) 



(13) 



(14) 



The constants Ci ( 2 are constants of integration. Since we 
know that the term linear in G does not change the field 
equations, we can without any loss of generality assume 
G\ = 0. Furthermore, if we wish to ensure that for to = 
2/[3(l + w)] and if = 4/[3(l + w) 2 ], the theory reduces to 
GR, i.e., f(G) = 0, then this constrains C2 = 0. Hence, 
for an exact power law solution to exist, the required 
form of the function / becomes 



f(G) = A m ^G + B, 



u gjrn(l+w) _ 



(15) 



We note that the above form of / identically satisfies the 
other field equations, if we similarly convert them as dif- 
ferential equations in G space. The coefficients A m , B mw 
are real- valued and non-zero unless m — 1 , in which case 
a(t) oc t. In general, the function f(G) is real- valued only 
if G [dim > 0, which is true by construction for the exact 
power law solution {8} . Similar results were found in [28| 
using a scalar-tensor approach to /(C?)-gravity. 

It is interesting to note that an exact GR-like solution 
(m = 2/3(l + u;)) is possible with non-zero C2 for f(G) — 

l+3w 



Co-existence of decelerating power— law solutions 
and accelerating solutions 

To mimic the standard expansion history of the uni- 
verse in /(C/)-gravity, we assume there exists an exact 
decelerating power-law solution, and that the universe 
was well described by this solution in the past, before 
coming to the accelerated phase. As argued in the previ- 
ous section, the existence of the exact solution fixes the 
form of f(G) as in (|15[) . and the deceleration fixes the 
power to as < to < 1, implying a m < 0. 

Now, if any additional accelerating solution exists in 
the whole solution space of the model, then G > for 
this solution as evident from (0 . However, we can see 
from the form of f(G), that this is not possible, as for 
G > the function is no longer real valued. 

This problem can be artificially remedied b y in clud- 
ing absolute values of G in (fT~5|) (similarly to [31(), i.e. 
redefining 



HQ) = K 



B„ 



~ jm(l+w) 

G 



(16) 



C2G 4 ( 1+W ) for values of w for which f(G) is real-valued. 



This function f(G) now seems to allow for both a de- 
celerating power-law solution and accelerating solutions. 
However, f(G) is not differentiable at G — and hence 
no longer a C 2 -function, which is required for any La- 
grangian in the Einstein Hilbert action. Any expansion 
history evolving from deceleration to acceleration must 
pass through G = 0, and we can easily see that the field 
equations are no longer defined at this point. Hence we 
can conclude that no well-defined C 2 action in f(G)~ 
gravity can allow for exact decelerating power-law solu- 
tions to co-exist with accelerating solutions. 



Comparison to f(R)— gravity 

It is interesting to note that in /(i?)-gravity the same 
argument cannot be carried through. One important dif- 
ference between f(R)- and /((J)-gravity is that in f(G)~ 
gravity the energy-momentum tensor decouples from the 
correction terms. In /(i?)-gravity on the other hand, the 
correction term modifies the matter energy-momentum 
tensor, which becomes rescaled to = jrr^j T^. Con- 
sequently it is not possible to convert the field equations 
in /(i?)-gravity to linear differential equations for the 
function f(R) in R space, but instead one obtains a 
non-linear differential equation. This equation cannot 
be solved easily, and in particular does not have a unique 
solution as is the case in the f(G) analogue. Additionally, 
since in /(Cy)-gravity matter decouples from the correc- 
tion terms any matter dominated era must strictly obey 
the corresponding GR type power-law behavior. 

The analogue also breaks down in the sense that any 
cosmologically viable model must make a transition from 
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deceleration to acceleration, and therefore go through 
Q = (see (J7J). This poses a problem for the simple mod- 
els including terms of the form Q n with n < 1, which are 
not differentiable at Q = 0. The value of the Ricci scalar 
R on the other hand does not necessarily pass through 
R = 0, and therefore models including terms like R n 
may still be viable. In fact if f(R) } B,\R^o± = for any 
/(i?)-gravity model, then the plane R = is actually an 
invariant sub-manifold. This means that no solution can 
cross that plane in phase space. For example, for any 
n < 1, R n is not differentiable at R = 0, but any solu- 
tion can only approach this plane from either side and 
not reach it in finite time. In fact, as shown in and 
[2fj| . in i?™-gravity a Friedmann like matter dominated 
decelerated solution co-exists with a de-Sitter like accel- 
erated solution. These two fixed points can be linked 
without R changing sign, i.e. the field equations are well 
defined and fully differentiable along the entire orbit. 

DISCUSSION AND CONCLUSION 



at the transition from deceleration to acceleration. 

One possible way around this problem stems from the 
fact that our conclusions are based on the requirement 
that there exists an exact power-law solution. If however 
the scale factor behaves like a(t) oc e At t m , then at early 
times a(t) oc t m . In this case the basic assumption of our 
analysis is not satisfied. 

We therefore conclude that when looking for cosmo- 
logically viable solutions or fixed points in the dynami- 
cal systems state space in /(^)-gravity, one should not 
look for exact power-law solutions, but rather exact solu- 
tions that approximate decelerating power-law solutions 
at early times. Furthermore, we must restrict ourselves 
to functions f(G) that are at least C 2 functions, and in 
particular differentiable at Q = 0. 

The authors would like to thank the National Research 
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The expansion history of the Universe is thought to 
have undergone a phase of decelerated power-law ex- 
pansion followed by late time acceleration. Therefore, 
power-law solutions play an important role in cosmol- 
ogy as matter dominated phases that later connect to an 
accelerating phase. 

In this work we have shown that exact power-law solu- 
tions in /(C/)-gravity only exist for the very special class 
of models given in (|15|) . In particular, many of the popu- 
lar /(^)-models, e.g. all examples considered in [32j |. do 
not allow for any exact power-law solutions. This means 
that for these models, there exists no exact matter dom- 
inated solution, and therefore these models may not be 
of cosmological interest. 

Furthermore, we have shown that decelerating power- 
law solutions cannot co-exist with any accelerating so- 
lutions, since no differentiable function f(G) can admit 
both decelerating power-law solutions and accelerating 
solutions. The problem may be remedied by includ- 
ing absolute values in (fT"5|) . However, this makes f(G) 
non-differentiable at Q = 0, which is the value Q takes 
when the scale factor evolves from deceleration to ac- 
celeration, and therefore cannot be ignored. This result 
seriously constrains the cosmological viability of f{G)- 
gravity, since it may not be possible to obtain an expan- 
sion history similar to the A - CDM model in this context. 

This issue has not been addressed in recent papers [3l| , 
where dynamical systems methods were used to study 
certain classes of /((?)-gravity models. It must be em- 
phasized that even if equilibrium points corresponding to 
both decelerating power-law and accelerating solutions 
are found to co-exist as in [3l|, there cannot exist any 
trajectories connecting these points, since these solutions 
can only co-exist for actions that are non-differentiable 
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